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Abstract 

This paper presents the measurement of the scattering of 172 MeV/c muons in 
assorted materials, including liquid hydrogen, motivated by the need to understand 
ionisation cooling for muon acceleration. 

Data are compared with predictions from the GEANT 4 simulation code and 
this simulation is used to deconvolute detector effects. The scattering distributions 
obtained are compared with the Moliere theory of multiple scattering and, in the 
case of liquid hydrogen, with ELMS. With the exception of ELMS, none of the 
models are found to provide a good description of the data. The results suggest that 
ionisation cooling will work better than would be predicted by GEANT 4.7.0p01. 
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1 Introduction 


The storage of intense muon beams is required for the construction of a Neu¬ 
trino Factory or muon collider. The ionisation cooling technique(P) is a leading 
contender to maximise the intensities for the Neutrino Factory and it is even 
more likely to be used at a muon collider. This technique relies upon the cool¬ 
ing effect of dE/dx losses in low Z materials coupled with R.F. acceleration 
to reduce the emittance of a beam. Multiple scattering in the material heats 
the beam and the overall cooling is thereby reduced. It is therefore important 
to conhrm the theoretical understanding and computational implementations 
of multiple scattering of muons of around 200 MeV/c. There is no published 
measurement of muon scatteriM in this region, and furthermore data on elec¬ 
tron scattering with similar (3^ suggest an excess in the tails compared with 
predictions from Moliere theory(j^, especially for low Z materials. This was 
first observed in ref. ( 0 ). 

Many problems involving multiple scattering have em^oyed GEANT 4 (j^ to 
evaluate the effects. A new model of scattering, ELMS(j^, has recently become 
available. This is a return to electromagnetic first principles, but often uses 
numerical integration rather than closed-form approximate integrals. 

The MuScat experiment was accepted by TRIUMF in 1999 to run in the Mil 
muon beam-line, and an engineering run took place in summer 2000. A new 
scintillating fibre tracker was built for the physics run. The data presented 
here were collected in April and May 2003, in the M20 muon beam-line in 
TRIUMF. 

Section 2 contains a description of the principle and technology of the MuScat 
experiment. This is followed in sections 3, 4 and 5 by details of the per¬ 
formance of the tracker, calorimeter and scintillator subsystems. The beam 
properties are described in section 6. These descriptions are accompanied by 
discussion of the GEANT 4 model which was used to deconvolute the detector 
effects, a process discussed in section 7. The data are presented in section 8, 
and conclusions can be found in section 9. 

The detector performance has been modelled using GEANT 4.7.0p01. This 
toolkit was used to make a complete description of the geometry shown in 
hgure 1. Gonsiderable attention has been attached to the simulation of the 
collimator system and the hbre tracker, which are critical to the experiment 
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Fig. 1. The MuScat detector, from the GEANT 4 model. 


and whose performances are sufficiently complex that they are hard to extract 
from the data alone. 


2 The Experimental Apparatus 


The technique of the MuScat experiment was to prepare a narrow collinear 
muon beam by collimation, and to use it to illuminate a variety of target 
materials. The position of the outgoing particles was measured after travelling 
a sufficient distance that the position and angle were highly correlated. 

The key components are the collimation system, a wheel with a selection of 
solid targets and an optional liquid hydrogen (LH 2 ) vessel, a scintillating hbre 
tracker, and the trigger scintillators. In addition to these the TINA (TRIUMF 
Iodide of Sodium) calorimeter (0) was useful to help understand the system. 
There were also some extra scintillators employed as veto counters. 

The passive components have been implemented in a GEANT 4 model in 
a simplihed manner, sufficient to simulate the materials which would affect 
the observed scattering distributions. The position of many externally visible 
components were checked with a precision of a few mm in the M20 beam-line 
in TRIUMF, and no deviations found. The position of the target wheel is 
known from the engineering drawings. The locations of the scintillators and 
the TINA calorimeter are known from the measurements. 
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The coordinate system chosen is that z is defined approximately by the beam 
direction, y is vertical and x is horizontal so as to define a right handed coordi¬ 
nate system. Since the collimator slits are narrow in the vertical direction, the 
y coordinate of a muon is the measure of its scattering. The precise definition 
of the z axis is that that the first two detector planes are both centred on 
X = y = 

2.1 The M20 beam-line 


The experiment used the M20 beam-line at TRIUMF. The extracted pro¬ 
ton beam from the cyclotron interacted in a target and produced pions at 
a production angle of 55 degrees. The beam-line had a muon decay channel 
consisting of two dipoles with quadrupole focusing. The beam-line was tuned 
to optimise the capture of high momentum muons from forward pion decays. 
The mean beam momentum was 172±2.0 MeV/c and the accepted momen¬ 
tum bite was approximately 1 MeV/c. The beam delivered to the experiment 
consisted of muons with a small pion background, but a mixed beam could be 
produced, as described later in the discussion of the time of flight system. The 
trigger rate was 80 particles/sec; the particle rate incident on the first scintil¬ 
lator, before collimation, was three orders of magnitude larger. The spot size 
at the first collimator was 50 mm in diameter. 


2.2 The Collimator 



Fig. 2. The second collimator piece, made of 160 mm of lead. On the left is the com¬ 
puter model, and on the right the two interlocking pieces of lead used to construct 
it. 


The collimator system consisted of two lead blocks with slits in them, aug¬ 
mented by other active and passive collimators. The accepted particles had to 
pass through both slits. This defined a beam which was narrow in one direc¬ 
tion, vertically, but much broader in the other to maintain rate. The upstream 
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lead block was 40 mm thick, with a simple slit, 20 mm by 2 mm. The second 
block, shown in hgure 2, had a trapezoidal slit, 12.8 mm by 2.0 mm at the 
upstream end, expanding to 20.0 mm by 2.88 mm at the downstream end. The 
collimator blocks were each made from two interlocking pieces of lead, with 
the lower one hxed to two steel rods by aluminium supports. The distance 
between the outer faces of the collimators was 963 mm. 

The system was behind a shield wall; 8 cm thick steel with a 4 by 3 cm opening 
aligned with the lead collimators. There were also eight auxiliary collimators, 
made from lead discs 20 mm thick and with 20 mm diameter openings, which 
were mounted on the steel rods between the two main collimators. They were 
intended to intercept any particle which might have scattered in the jaws of the 
hrst collimator and then missed the second one entirely. They were designed 
so particles scattered from them could not pass through the slit in the hnal 
collimator without another interaction. The entire tube was wrapped in lead 
sheet to a thickness of 6 mm. This ensures that all muons reaching the tracker 
have passed through the entire collimator system. 

Finally, the downstream passive lead collimator was supplemented with two 
scintillators, 3 mm by 3 mm by 30 mm. These were positioned just at the 
front face of the collimator, above and below the slit. They were designed to 
record particles which might otherwise have hit the lead and then emerged 
into the slit with a large scattering and energy loss. They were read-out by 7 
of the clear 1 mm hbres used for the Sci-Fi tracker. The hbres were attached 
to each end and taken to a photomultiplier tube (PMT). 

The width of the collimated beam must be well simulated, and unfortunately 
this is not the case with the assumption of a perfect geometry. The apparent 
slit width and also the observed tails are both wider in data than in the 
simulation. The only credible way to increase the collimated beam width is 
to hypothesise that the second collimator was not perfectly closed. It was 
constructed as two blocks of lead which dovetail neatly and hrmly but are not 
attached together. The lower block was held hrmly to the rails, and the whole 
assembly is prevented from rotating by a grub screw through the collimator 
tube attached to the top lead block. If the system was somehow pressured after 
or during installation the grub screw could pull the top lead back, widening the 
jaw. The simulated distribution of particles through the expanded collimator 
assumption is shown in hgure 3, and this is used hereafter. 

This is unfortunate, as not only does the central distribution widen but the 
inner faces of the collimator can become directly exposed to incident particles, 
allowing scattering from them. These tails are indeed observed. The energy 
distribution of the tail particles, as measured in TINA, is consistent with what 
would be expected if this is the correct explanation. The size of the movement 
is small. The best representation of the observed data is obtained if the jaw is 
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Fig. 3. The distribution of particles at detector plane 1 in data and simulation for 
no target. The error bars include systematic error estimates. The large error bar 
just to the right of the peak includes the effect of possible misalignment. 

opened from 2.0 mm to 2.48 mm. 


2.3 The target system 


There were two separate target systems: a wheel with a choice of solid targets 
and a vessel for liquid hydrogen. The target wheel was controlled by a stepper 
motor, and presented a choice of twelve target positions to the beam. The 
position was monitored by an optical LED system which allowed readout of 
the target number and a confirmation that the position was exactly correct. 
Table 1 shows the targets used, and the number of events selected for analysis 
from each. The density and radiation length, Xq, are also shown. 

The LH 2 vessel was crucial to the experiment, because scattering in liquid 
hydrogen was the most important goal. This was designed and built by the 
cryogenic targets group at TRIUMF. There is a contribution to the scatter¬ 
ing from the containing windows, and therefore it it is essential to minimise 
and measure this accurately. The vessel could be filled with liquid or gaseous 
hydrogen, and the measured distributions with gaseous hydrogen are used to 
characterise the windows. The temperature remained between 15 and 23 K 
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No. 

Material 

Thickness, mm 

Density, g/cc 

Xo, % 

Events 

0 

Lithium 

12.78 

0.53 

0.82 

805239 

1 

Lithium 

6.43 

0.53 

0.41 

1224756 

2 

Lithium 

6.40 

0.53 

0.41 

882449 

3 

Lithium 

12.72 

0.53 

0.81 

1215336 

4 

Beryllium 

0.98 

1.85 

0.28 

500766 

5 

Beryllium 

3.73 

1.85 

1.06 

1186528 

6 

CHs 

4.74 

0.93 

0.99 

802426 

7 

Carbon 

2.50 

1.69 

1.53 

801899 

8 

Aluminium 

1.50 

2.70 

1.69 

1201280 

9 

None 

- 

- 

- 

2259476 

10 

Iron 

0.24 

7.86 

0.82 

845020 

11 

Iron 

5.05 

7.86 

28.68 

1225435 


Empty H 2 

109.0 

- 

- 

2694511 


H 2 

109.0 

0.0755 

1.31 

2267683 


Empty H 2 

159.0 

- 

- 

812730 


H 2 

159.0 

0.0755 

1.90 

1127045 


Table 1 

The targets, and the number of data events selected for each, 
for these cycles, and there was always a slight overpressure. 

The vessel was 150 mm by 100 mm by 100 mm. The 150 mm length was 
bored with an 8 cm diameter hole, and an identical hole was made through 
from one side. These were capped by four identical windows in the sides and 
ends, so that by rotating it a choice of hydrogen thicknesses was presented. 
The windows bulged by 4.5 mm, with an estimated error of 0.3 mm, making 
the total length 109.0±0.6 mm or 159.0±0.6 mm. They were made from mylar, 
5 thousands of an inch thick, wrapped in 6 layers of super-insulation. 


3 The Sci-Fi tracker 


The scintillating hbre (Sci-Fi) tracker comprised three double hbre planes, 
read out by twenty-four HPK multi-anode PMTs (@). The planes were about 
1100 mm from the target. The hrst and second planes were mounted together, 
which means the active elements were displaced by 50 mm. While for the third 
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plane the displacement was 66.9 mm. 



Fig. 4. The layout of the scintillating fibres. Each plane consists of both horizontal 
and vertical fibres, overlapping in two rows to give complete coverage. In each row 
the distance between fibre centres was 1.17 mm and the distance between the rows 
was 1.17 mm. They are potted in black epoxy. 

The donble planes consisted of two planes (each constrncted as shown in £g- 
nre 4), one comprising 512 horizontal hbres and the other the same nnmber of 
vertical hbres. They have a pitch of 0.585 mm, covering 300 mm^. The hbre 
diameter was 1 mm, inclnding the onter 4% cladding. They were arranged in 
donble rows, as in hgnre 4. They were held in a Delrin frame, 350 mm sqnare 
with a 300 mm sqnare central area. 2048 holes were drilled with a CNC ma¬ 
chine to dehne the locations of the two ends of each hbre. The hbres were held 
nnder tension and potted in black epoxy. Finally the ends of the hbres were 
cnt. 



Fig. 5. The physical arrangement of the Sci-Fi tracker. The left hand figure shows 
one plane mounted inside the vessel. The black frame at the top encloses the scintil¬ 
lating fibres, and the clear routing fibre bundle can be seen. The right hand figure 
shows the front faces of the fibre bundles, with a PMT mounted onto one of them. 

Each double plane had four Delrin edge connectors, made in a similar manner, 
and drilled in one operation with the frame so that the holes matched. This 
can be seen in hgnre 5 (left). Into each of these connectors were potted 512 
double-clad 1 mm radius clear hbres which were to guide the light to photo¬ 
multipliers. This was done in two groups of 256, and the far ends of the hbres, 
after about 300 mm, were arranged into a square 16 by 16. These were potted 
in an aluminium former which was to also act as a vacuum seal for the detector 
vessel. The natural tendency of the hbres to close-pack made this construction 












difficult, and the edge fibres were not always perfectly placed. The best bundles 
were selected for the first detector and for measuring y. The 16 mm square 
fibre bundles presented on the outside of the pressure vessel were connected to 
the 16-anode PMTs, as shown in figure 5 (right). The readout was done with 
a custom-made sample and hold system, which provided a stable, low-noise 
readout for the system. 

Each set of 256 fibres is read out at both ends. One end of the fibres is read 
out by a PMT which measures the signals from 16 adjacent fibres in each of 
its (16) anodes. The other set of ends is distributed across the PMT surface 
so that each PMT anode reads every sixteenth fibre. In this way, one PMT 
is used to measure which group of 16 fibres has a signal and the other which 
fibre within the group. This multiplexing scheme can only work with a single 
charged particle present. 

Simulation of the tracker was fairly complete. The light yield per fibre in the 
simulation was tuned to the data using the mean measured signal in identified 
hits using the thick iron target. For a few percent of fibres with low light yields 
this does not converge satisfactorily; for these the rate of hits is used instead. 
Cross-talk and signal shape were simulated due to the light spreading in the 
1.5 mm glass at the front of the PMT, the photo-electrons jumping to other 
dynodes, a leaky dynode chain where an electron can skip a stage or move 
into the neighbouring chain and about one percent of electronic cross-talk. 




Fig. 6. Left: The number of ADC counts on every channel in data and simulation 
with the thick iron target. The peak at 600 is due to single photo-electrons. Right: 
The distribution of the number of fibres with a signal of over 0.2 photo-electrons 
observed in data and simulation. 


The distribution of ADC counts recorded on all channels is in figure 6 (left). 
The pedestal is clearly visible, with a peak due to single photo-electrons at 
around 600 counts which was used to measure the gain. The upper edge is cut 
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off by the saturation of the ADCs. The simulation includes these effects and 
shows reasonable agreement. 


The distribution of the number of reconstructed hbres with signals in data 
and simulation is shown in hgure 6 (right). The data has events with no 
hbres hit which are not reproduced by simulation. These are trigger accidental 
coincidences. There are also events with several times the normal number, 
which arise from multiple tracks. Both types of events are removed in analysis. 
The central region has a reasonable representation. 
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Fig. 7. The distribution of particles at detector plane 1 in data and simulation 
when the thick iron target is in place. 


The recognition of a hit in the tracker depends upon the observation of at 
least one photon at either end, and as the average light yield is around four 
photons then any poor hbres with lower yields will have signihcant inefficiency. 
This is examined in hgure 7, which shows the results for a thick iron target, of 
28% Xq. The scattering distributions in this target are dominated by the cen¬ 
tral Gaussian core, and the simulations are therefore expected to be reliable. 
The purpose of this target is to illuminate the entire detector with tracks for 
use in alignment and efficiency studies. There is a distinct drop in efficiency at 
the very centre of the distribution, and also some structure at about -10 mm. 
The distribution is shown for the entire detector, and the agreement is well 
within 10% in most bins, with the centre and the edge bins showing differences 
around 20%. It seems that the simulation of the edges of the PMTs is some¬ 
what defective, probably due to the difficulty of locating all the clear hbres in 
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the correct places at the front of the PMT. 


4 The TINA Calorimeter 


An existing calorimeter from the TRIUMF Laboratory, ‘TINA’(0), was used to 
understand the beam properties. This is a cylindrical sodium iodide calorime¬ 
ter, 460 mm diameter and 510 mm long, split into two crystals in depth. It 
is read out using 7 PMTs on the back-face, one in the centre and the other 
six in a circular pattern around it. The calorimeter is contained in a metal 
blockhouse to reduce external backgrounds. This blockhouse has a small lip 
that overlaps the outer edge of the Nal by about 40 mm. The device has a 
measured energy resolution (fwhm) of 3.6% at 90 MeV for electrons, with an 
energy dependence of 

As TINA was too small to cover the complete angular acceptance of MuScat, 
it was mounted with the centre shifted upwards by 70 mm with respect to 
the axis of the experiment. This allowed a complete angular coverage in the 
positive y direction. In addition, the energy deposition was integrated over 
two time periods: tina_early, corresponding to 500 ns from the SIB trigger, 
and tina_total, corresponding to 10 /iS from the trigger. In general, the hrst 
of these will have only the kinetic energy of the muon, while the second also 
contains the electron from the decay of the muon, which has come to rest in 
TINA 


30000 
25000 
20000 
15000 
10000 
5000 
0 

200 300 400 500 200 300 400 500 

Fig. 8. The tina_early ADC signal in the TINA calorimeter for events selected by 
timing cuts (see section 6.2) to be either mixed pion and muon events (left) or the 
muon sample used for analysis (right). A requirement on tina_total - tina_early 
has been made so that the muon decay is known to occur after the tina_early gate 
closed, but a small decay tail can still be seen. The peak at 245 ADC counts is the 
pedestal. 

The performance of TINA is illustrated by hgure 8. The pedestal is at 245 
ADC counts. The left hand plot (pion sideband), shows the separation between 
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pions and muon peaks. Decay electrons have been supressed as discussed in 
the caption. If this had not been done there would have been larger tails, and 
an additional 0.4% would have been assigned to the background. The TINA 
calorimeter was modelled in GEANT 4, but these simulations are not used in 
this analysis. 


5 The Scintillators 


There were seven instrumented plastic scintillators used in the experiment. 
Three of these made up the trigger system, and the other four were veto 
counters. 

The trigger was provided by a logical and of two overlapped huger scintilla¬ 
tors, 3 mm square in area and 30 mm long, known as SIA and SIB. These 
were mounted in front of the collimator to record particles before entering the 
system. They gave a precise timing signal, used for the time of hight mea¬ 
surements discussed later. This was taken in coincidence with the signal from 
a 330 mm square scintillator, referred to as S2, which was read out by two 
PMTs, top and bottom, whose results were combined as a logical or. S2 was 
placed after the Sci-Fi tracker and before TINA as shown in hgure 1. 

Two large area veto scintillators were butted together just in front of SI; 
one had a notch removed that corresponded approximately to SI. These are 
referred to as VI. The other two veto scintillators were bars 3 mm by 3 mm 
by 30 mm which were mounted on the front face of the second lead collimator, 
just above and just below the slot. 

The common start for all the TDCs (Time to Digital Converters) in the ex¬ 
periment came from the SIB scintillator. A delayed version of SIB was then 
read by a TDC. This should have given a hxed time, but some events were 
several ns different and they were removed by offline selection. An additional 
signal from the cyclotron, repeated each 43 ns period, provided a useful timing 
reference. 

The GEANT model included accurate descriptions of the geometry of the scin¬ 
tillator system. Electronic handling of the signals, in terms of pulse shaping, 
and discriminator threshold for the TDCs was included in an approximate 
way so that time-walk effects were reproduced. 

In order to reproduce the data, the huger scintillators comprising the active 
collimators had to be positioned so that they overhung into the collimator gap 
slightly: the top by 190 /xm and the bottom by 70 //m. This reinforces the idea 
that the collimator slot had somehow opened. 
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6 Beam properties 


There has been no attempt to simulate the beam-line. However, the x dis¬ 
tribution of the observed beam has a very small dependence upon the initial 
beam conhguration, and for y it is totally negligible. The beam divergence is 
instead determined by multiple scattering in the SIA and SIB scintillators, 
which guarantees a minimum spread, followed by the collimation system which 
limits the maximum angle. 


6.1 Beam momentum determination 



200 300 400 500 600 700 800 900 1000 


TDC counts 

Fig. 9. The time of arrival of the beam cycle signal relative to the SIB scintillator 
for run 177, which was tuned to have a mixed particle content. The four lines 
guide the eye to pions, muons, protons and electrons. The open histogram is for all 
events that are in time on SIB and the hashed area for events that meet a tighter 
requirement on SIB within 0.25 ns of expected, which excludes protons. A double 
peak structure can be seen for the precise peaks; this is not understood. 

A special run was taken with the momentum of the second dipole in the 
beam-line turned down by 2%. This meant that the forward muons normally 
accepted have a large contamination from other species. The presence of multi¬ 
ple particle types has been used to measure the beam momentum. The time of 
arrival of the cyclotron signal relative to SIB is shown in hgure 9. The major 
peaks are at 500 (pions), 650 (muons), 750 (protons) and 900 (electrons). In 
general the fastest particles give the largest times, because what is measured 
is the arrival time of the cyclotron signal after SIB. In the case of protons, the 
velocity is only 0.2 c, and they arrive seven beam crossing periods late. The 
width in the electron and pion signals are similar, but the others are broader. 
For protons this is because the very low velocity leads to a much larger spread 
in time, while the muon arrival time depends upon the decay point. There is a 
distinct peak of muons which were produced before the hrst bending magnet. 
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Mass, MeV/c^ Mass, MeV/c^ 

Fig. 10. Fitting the beam momentum from the observed peaks using the cyclotron 
signal. The negative of the arrival time of the cyclotron signal is plotted; for protons 
a point is shown both at the measured time and delayed by seven beam cycle periods. 

The momentum of the beam has been estimated by htting the arrival times 
of all four particle species. Three parameters are required for this: the time 
offset, (dehned by the arrival time of a particle going at the speed of light), the 
length of the beam-line and the beam momentum. The results are shown in 
hgure 10, where a momentum of 168.9 MeV/c is found and the four points are 
seen to be consistent. This run was taken with the hnal dipole held reduced 
by 2%, so under normal running conditions the beam momentum is 172±2 
MeV/c. 

The spread in momentum can be estimated by noting that it affects the veloc¬ 
ity of pions but not electrons. The quadrature difference in width of the pion 
and electron peaks calculated as in hgure 9 but after accounting for small 
time-walk ehects gives the spread in momentum as being between 1.0 and 

1.2 MeV/c. 

6.2 Beam composition 


The recorded sample under normal conditions has a large proportion of muons. 
This can be seen in hgure 11. Protons are ranged out in the SIA trigger, and 
not generally recorded. This hgure shows a somewhat worse time resolution 
than in hgure 9, which is not understood but was constant for the data used 
for the main analysis. 

The muon fraction is further raised by demanding that the recorded time of the 
TOF signal is within the section labelled pi in the hgure. The energy recorded 
in TINA for the subset of events within its acceptance is used to evaluate the 
purity of this sample. The left hand plot on hgure 8 (pion sideband), has two 
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TDC counts 


Fig. 11. The time of arrival of the beam cycle signal relative to the SIB scintillator 
for run 290, an ordinary run with no target. The graph corresponds to events where 
the digitised time from SIB was correct. The section labelled fi contains those events 
used for the scattering analysis. The small peak to the right contains electrons. The 
resolution is slightly worse than in figure 9; this is not understood. 

large peaks due to pious (left) and muons (right). The distribution is htted 
using two Gaussian distributions to determine the peak positions and widths. 
The relative amplitudes of these contributions are then htted to the right-hand 
plot, where a 0.8% pion component is found. Some contribution to the second 
Gaussian comes from resolution tails, and therefore the pion contamination is 
almost certainly overestimated. 


7 Analysis Method 

7.1 Event Selection 


The selection of events for analysis requires that they pass a number of se¬ 
lection criteria. Firstly the high voltage must have been nominal and the 
data-acquisition in a normal state. Then we require that the event was not 
taken during a period identihed as having a problem with the time of hight 
system. 

Next the number of reconstructed points on the hrst plane of the tracker is 
required to be one or two; more hits would make the position of the muon very 
ambiguous. Then the summed signal from the two internal veto scintillators 
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had to be consistent with the pedestal, and the time from SIB within a few 
counts (0.2 ns) of nominal. 


The hit with the largest light yield is taken as the estimated particle position 
at each detector plane. Finally, the positions of the hits on the hrst and second 
planes are required to be consistent with the hypothesis of a particle coming 
from the target. The second plane is about 4% further from the target, so it 
is required that 0.96 times the y position on the second plane be within 5 mm 
of the y position on the hrst plane. This procedure is simpler than performing 
a full track £t, and less dependent upon alignment systematics. 


1.2 Deconvolution technique 


It is possible to compare the distributions of hits at the hrst detector plane. 
However, the more fundamental quantity is the angle through which the parti¬ 
cle scattered in the target, and the ehects of beam width, efficiency, resolution 
etc. must be corrected for to see this. This deconvolution is done using the 
GEANT 4 simulation, with the following approximate formula: 

D = B + Dt, + R-€-Q ( 1 ) 


Where D is the observed position data, B a background of particles not passing 
through the target, is the contamination from pions, R is the response of 
the detector to a particle dehected through angle 6y, e is the efficiency of 
the detector for particles dehected through angle 6*^, and 0 is the projected 
scattering distribution in the target. 

The background, H, is found from simulation as those muons which managed 
to meet the trigger conditions without passing through the target. There are 
typically 0.125% background events. The pion contamination, is taken 
from the data pion sideband. The default value is 0.8% pions, and as the 
sideband is approximately 50% pions, a 1.6% admixture of the data sideband 
is added to the simulation. 

The response and efficiency matrices, R and e are taken from simulation. The 
efficiency is found by running GEANT 4 twice. The hrst run is a simulation of 
the target only, illuminated by a monochromatic collinear beam of muons. This 
is used to hnd the true 0 distribution in GEANT 4. The second run is the full 
simulation, including trigger and tracking acceptance cuts. The distribution 
0 for those muons accepted is calculated, and the ratio of these is e. In the 
above formula it is treated as a diagonal matrix. Finally, R is hlled as a 2D 
matrix, with each accepted simulated event being entered at a point given by 
its dy and the y position at the hrst measurement plane. It therefore gives 
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the probability of a given measured position for each possible true deflection 
angle. 

The equation is solved for 0, imposing a requirement of symmetry about 6 = 0, 
using Minuit(j^, which also hnds errors and correlations. It has been checked 
that the technique is mathematically correct, in that if given simulated data 
as an input the scattering distribution from that simulation is recovered. 

This formalism does not explicitly include the effect of energy losses and scat¬ 
tering in the X direction. The matrices R and e therefore depend on the 
simulation of multiple scattering and energy loss in the GEANT model. In 
particular, the probability that a particle completely misses the detector in 
the X direction as a function of its Y angle must be taken from simulation. 

The problem of instabilities in the deconvolution has been dealt with by reduc¬ 
ing the number of bins in the deconvoluted data, where 21 are used, compared 
to the raw data, which has 57. Furthermore the condition of symmetry is en¬ 
forced in the £t, so there are essentially only 11 bins. Finally, the outermost 
bin in true scattering angle runs from 115 mrad to tt, but is very weakly con¬ 
strained by the data as the majority of events in this bin are not recorded. 
It also has a strong dependence on the simulated prediction of the fraction 
recorded. It is not shown in the plots which follow. 


8 The results 

8.1 Comparisons at detector level 


Data taken with no target, an empty hydrogen vessel or a very thick iron 
target are not of great intrinsic interest, but were necessary to understand the 
detector. The hrst, the bare collimator distribution, is shown in figure 3. The 
hnal results are based upon deconvoluting the data using the simulation, and 
so discrepancies here will translate into systematic errors on the results. In 
the case of hydrogen data, the equivalent of the bare collimator is the empty 
target vessel. The distributions obtained for this in data and simulation are 
showing in hgure 12 for the case of the longer target conhguration. 

The error bars shown in these plots are statistical only for the data but are 
a combined statistical and systematic error on the simulation points. This 
systematic error explicitly includes the difference between these plots, and so 
by construction all ratios must be compatible with one to within one combined 
error. 
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Fig. 12. The distribution of particles at detector plane 1 in data 
the empty H 2 vessel in the 150 mm position. 


Position, mm 

a and simulation 


for 


The other control plot is figure 7, which shows the results for a thick iron 
target, of 28% Xq. The agreement is well within 10% in most bins, with the 
centre and the edge bins showing differences around 20%. As the edge bins 
have also got a large background subtraction, the region of the detector used 
for analysis and subsequent hgures is restricted to the central 270 mm. 


The raw distribution at detector level for a few of the targets are now shown. 
Figure 13 shows the combined data from the two thin lithium targets. The 
thickness of this sample corresponds to the second lowest fraction of a radiation 
length of material used, and it therefore has a narrow distribution. GEANT 
tends to overestimate the distribution, starting from about 10 mrad. Note 
that outside about 100 mm the distribution is completely dominated by the 
scattering from the collimator, see hgure 3, and the contribution from the 
lithium is negligible. Thus the convergence of data and simulation in the tails is 
more related to the simulation of the collimator than the physics of scattering. 

The thin iron target, hgure 14, is the highest Z target studied. The GEANT 
description is better than for any of the lighter targets, but there is here a 
small region below 40 mm where the scattering distribution is underestimated 
which is seen in no other target. 


The hnal comparison presented at the level of the detector is for the longer 
hydrogen target. This resembles the lithium target, but the differences between 
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Fig. 13. The distribution of particles at detector plane 1 in data and simulation 
for thin lithium, combining both targets. This distribution is used to check the 
agreement of our data with GEANT. 


data and simulation are more pronounced, with the observed distribution less 
than half the predicted level at ±40 mm. 


8.2 Deconvoluted results 


The deconvolution procedure described in section 7 is applied to each of the 
targets. A series of systematic effects are allowed for, and plots for many of 
the targets are presented in hgures 16 to 22. In each case the statistical error is 
indicated with an inner error bar, and the combined statistical plus systematic 
error is indicated with the outer error bar. In most cases the statistical errors 
are only important for regions of the distribution where the probability of 
arrival is below about 0.1 per radian. The data for all targets are shown in 
table 2. 

The plots have two smooth curves superimposed: the solid line represents the 
Moliere calculation with a Z x (Z ± 1) factor for the large angle scatters, as 
proposed by Bethe (fli^ . and the dashed line represents the same calculation 
using a term. The difference corresponds to scattering solely on the nucleus, 
or also on the electrons. Note that the kinematics of scattering prevent a muon 
from making a single scatter from a free electron through an angle greater 
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Fig. 14. The distribution of particles at detector plane 1 in data and simulation for 
thin Fe, target 10 
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Fig. 15. The distribution of particles at detector plane 1 in data 
159 mm of liquid H 2 . 
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than me/rrif^, or 5 mrad. We therefore expect that for the description of the 
tail region, above 15 or 20 mrad, the term will be more appropriate. 



Fig. 16. The projected scattering angle distribution in data and simulation for thin 
iron, target 10. 

I Thin Iron. The data are shown in hgure 16. The two Moliere models 
are almost indistinguishable, and provide a good description of the 
data. In general, the GEANT 4.7.0pl simulations are somewhat in 
excess of the data above 35 mrad. Errors are around 5%, growing to 
30% in the hnal bin shown. 



y angle, radians 

Fig. 17. The projected scattering angle distribution in data and simulation for 
aluminium, target 8. 
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II Aluminium. The data, figure 17, are rather similar to the iron, with 
an overestimation of the tail by GEANT 4 of around a factor two. 



Fig. 18. The projected scattering angle distribution in data and simulation for 
carbon, target 7. 

III Carbon. The data, shown in hgure 18, reveal a visually apparent 
distinction between the Moliere models. Furthermore, the GEANT 
simulation appears to have an excess compared with data starting 
at lower angles, around 20 mrad, than it did for iron. 

IV Polyethylene The data for this target resemble beryllium, with a 
similar overestimate from GEANT of the scattering. They are not 
shown. The largest systematic contributions come from the simu¬ 
lated collimator description. 

V Beryllium The thick target data are shown in hgure 19. The GEANT 
model overstates the tails by approximately a factor of two. Of the 
two Moliere models, the is slightly preferred. No distinction was 
possible for materials with higher Z as the models were similar. The 
thin beryllium has larger errors, owing to the large background sub¬ 
traction required there, and is not displayed. 

VI Lithium There were two pairs of lithium targets with very similar 
thicknesses. The data for the thicker pair are combined in hgure 20. 

The Moliere model with a Z^ term provides a reasonable description 
of the data, but the Z x (Z -h 1) Moliere model is now clearly dis¬ 
favoured. GEANT continues to have tails about a factor two above 
the data. The data from the thin lithium targets have large system¬ 
atic and statistical errors in the tails, owing to the large fraction of 
background which has been accounted for in the deconvolution, and 
are not shown. 

VII Hydrogen For 109 mm of hydrogen it can be clearly seen in hg- 
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Fig. 19. The projected scattering angle distribution in data and simulation for thick 
beryllium, target 5. 



Fig. 20. The projected scattering angle distribution in data and simulation for thick 
lithium, both targets combined. 

ure 21 that neither Moliere model accurately describes the data. The 
model understates the region below around 25 mrad, and above 
60 mrad, although here the errors are much larger. GEANT 4.7.0p01 
predicts too much scattering, by a factor of around 4 at 50 mrad. 

The data from the ELMS model (j^ are also compared with our 
results. These seem to show good agreement at all but the high angle. 

This region is has the largest systematic errors. 

It can be clearly seen in hgure 22 that neither Moliere model accu- 


23 











































































Fig. 21. The projected scattering angle distribution in data and simulation for 
109 mm of liquid H 2 . 



Fig. 22. The projected scattering angle distribution in data and simulation for 
159 mm of liquid H 2 . 

lately describes the data for the long hydrogen target. The model 
seems correct beyond around 40 mrad, but in the multiple-scattering 
dominated region the Z x (Z -h 1) model is much more appropriate. 
GEANT 4.7.0p01 predicts too much scattering, but close examina¬ 
tion reveals that the excess is less than for the thinner sample. 

Overall, we see that when the two Moliere models are distinguishable, i.e. when 
Z is small, the data is described by the Z x (Z + 1) in the Gaussian core, while 
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the tails show better agreement with the version. This emphasises that the 
difficulty of describing the data is greater at low Z. 

GEANT 4.7.0p01 typically predicts about a factor two more scattering tail 
than is observed. This is not a very large effect given that the data cover 4 
orders of magnitude. However, in the case of liquid hydrogen the discrepancy 
grows to be a factor 4. 

The ELMS simulation describes hydrogen very well. 


8.3 Systematic Errors 


Systematic errors are evaluated by performing the complete analysis with 
different treatment of the data. The difference with the main results is taken 
as an estimate of systematic error. All errors are combined in quadrature. 


8.3.1 Collimator Description 

Three different techniques are used to evaluate the impact of differences be¬ 
tween the collimated distributions seen in hgure 3. The largest systematic 
errors are generally from uncertainty in our understanding of the collimated 
beam. 

A-I Collimator difference The bare collimator distribution seen in 
simulation is subtracted from the observed data. This residual is 
then added to the simulated distribution for whichever target is cur¬ 
rently under study. This is not completely correct, as it ignores the 
scattering in the target, but it provides an estimate of the systematic 
effects. 

A-II Collimator deconvolution The deconvolution to hnd the under¬ 
lying scattering distribution is performed on the data with no target 
present. This would yield a delta function, but does not due to bin 
width effects and discrepancies between data and simulation. The 
same is done in simulation, and the difference of the two is found. 

This difference is due to poor modelling of the collimated beam. 

The deconvolution for any given target is then repeated with this 
difference distribution included as an extra convolution term on the 
right hand side of equation 1. Any changes in the results are taken 
as systematic errors. 

A-III Internal veto scintillators The collimation system, described in 
section 2.2, has two internal veto scintillators mounted just before 
the second collimator. The data show a good separation between 
pedestal and hits, but there is no complete cross-check on this. The 
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simulated efficiency is lowered to 90%, (which also provides a slightly 
better description of the scattering distributions) and a systematic 
error is generated. 


8.3.2 Detector alignment 

Detector misalignment relative to the simulation will cause errors in the de- 

convolution. 

B-I Detector X alignment Any credible X misalignment will have 
negligible impact. 

B-II Detector Y alignment A transverse shift of the detector of 0.1 mm 
is allowed as a systematic error. The size comes from inspection of 
hgure 3. 

B-III Detector Z alignment The Z position of the detector affects the 
widths of the distribution. From a combination of studying internal 
distributions in the tracker and knowledge of the detector construc¬ 
tion a systematic error of 5 mm has been assigned. 


8.3.3 Detector modelling 

C-I Thick target comparison The muon tracking efficiency maps di¬ 
rectly into the scattering distribution, and as such it must be well 
controlled. The distributions in data and simulation observed with 
the thick iron target present, hgure 7, are used to monitor this. It is 
conservatively assumed that all differences are due to efficiency, and 
so the ratio of these distributions is used as an efficiency correction. 
The change in the unfolded results when this is applied is regarded 
as a systematic error. 

C-II Cross-talk suppression The default cut selection allows one or 
two hits to be reconstructed on the hrst detector plane. When we re¬ 
quire exactly one the probability of picking the wrong one is reduced. 
The analysis is repeated with this requirement and the difference is 
taken as a systematic error. 

8.4 Pion contamination 


As discussed in section 6.2, a pion background of 0.8% is estimated. However, 
this comes from the observation that hgure 8 has tails to the left of the muon 
peak; some of this background will arise from sources other than pions. The 
most extreme alternative assumption is that there are no pions present, and 
this is used to estimate a systematic error. 
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8.5 Deconvolution matrix 


The matrix used to relate true scattering angle to measured position is taken 
from simulation. It depends upon the physics in that simulation, especially 
scattering in the X direction. This dependence is reduced by re-weighting all 
the simulated events as a function of the space angle scatter with an ad- 
hoc simple function which moves GEANT 4 towards Moliere scattering 
above about 40 mrad. The central values of all the fits have been obtained 
using this procedure. A systematic error is assigned by comparing with the 
unweighted GEANT, except in the case of hydrogen where we have the ELMS 
simulation, and the alternative distribution is GEANT 4 re-weighted to match 
the predictions of that. This is typically the largest systematic error at the 
edge of the deconvoluted distributions. 


9 Conclusions 


These data allow the validation of codes to calculate multiple scattering in the 
regime relevant for ionisation cooling. In general Moliere scattering provides a 
good description, but for Z below about 4 there are differences which may be 
important. In particular, for hydrogen neither Moliere model used is appro¬ 
priate. The ELMS simulation provides a good prediction of the deconvoluted 
hydrogen distributions. 

The data show noticeable deviations from the predictions of the GEANT 4.7.0p01 
simulation code, which tends to overstate the scattering tail by about a factor 
of two, and more strongly the lower the Z of the material. For hydrogen it is 
about a factor of 4. 

The discrepancies observed previously (j3) are not confirmed. Indeed, the scat¬ 
tering predicted by the GEANT simulation needs to be reduced to match the 
data, not increased. 
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Xq Upper edge of bin, radians 

% Type 0.00269 0.00895 0.0162 0.0248 0.0347 0.0463 



data 

89.4+4.4 

37.3+1.7 

3.49+0.27 

0.31+0.05 

0.41 


G4 

75.6+0.1 

39.7+.04 

5.64+0.01 

0.56+.004 


data 

61.7+3.4 

39.3+1.3 

10.3+0.21 

1.14+0.06 

0.82 


G4 

53.1+.04 

38.3+.03 

12.9+0.02 

2.13+.007 


0.08±0.02 0.017±0.009 

0.13±.002 0.046±0.001 

0.19±0.02 0.053±0.009 

0.36±.003 O.lOliO.OOl 



data 

112.+4.9 

29.9+2.0 

1.43+0.17 

0.16+0.04 

0.28 


G4 

95.1.+0.1 

35.4+.03 

2.55+0.01 

0.33+.003 


Data 

54.2+5.4 

38.4+1.2 

12.7+0.63 

1.82+0.17 

1.06 


G4 

48.0+0.05 

36.5+0.04 

14.6+0.02 

3.10+0.01 


0.04±0.014 0.013±0.011 

0.09±0.002 0.033±0.001 

0.27±0.025 0.069±0.010 

0.60±0.004 0.170±0.002 


CH2 0.99 


data 

G4 


55.4±3.1 

50.2±.05 


38.2±1.3 

37.2±.04 


12.6±0.18 

13.9±0.02 


1.79±0.10 

2.68±.009 


0.28±0.020 0.073±0.010 

0.49±0.004 0.140±0.002 


data 55.7±3.0 37.8±1.3 12.7±0.18 1.97±0.08 0.29±0.022 0.076±0.010 

G4 50.3±.05 37.2±.04 13.7±0.02 2.67±.009 0.52±0.004 0.153±0.002 


A1 


data 41.7±2.4 33.Oil.0 17.2±0.20 4.82±0.10 0.96±0.040 0.198±0.014 

G4 39.8±.04 32.9±.03 17.0±0.02 5.14±0.01 1.15±0.006 0.327±0.003 


Fe 0.82 


data 

G4 


45.4±2.8 

44.4±.05 


34.4±1.1 

35.3±.04 


16.3±0.15 

15.8±0.02 


4.19±0.10 

3.64±0.01 


0.76±0.034 0.171±0.012 

0.78±0.004 0.250±0.002 



data 

49.5+2.7 

35.8+1.5 

15.7+0.41 

2.70+0.39 

1.31 

G4 

38.7+.04 

32.2+.03 

17.4+0.02 

5.33+0.01 


ELMS 

49.6+.07 

37.7+.05 

14.5+0.03 

2.43+0.01 


data 

40.6+1.7 

32.3+0.9 

18.3+0.18 

5.10+0.38 

1.90 

G4 

31.4+.04 

27.9+.03 

18.3+.025 

8.06+0.02 


ELMS 

40.9+.06 

33.8+.05 

17.5+0.03 

4.86+0.02 


0.34±0.061 

1.00±0.005 

0.26±0.004 

0.88±0.014 

2.22±0.008 

0.72±0.006 


0.061±0.012 

0.205±0.002 

0.052±0.001 

0.117±0.028 

0.466±0.003 

0.103±0.002 


Xq Upper edge of bin, radians 

% Type 0.0597 0.0754 0.0938 0.1151 3.141 


0.41 

0.82 

data 

G4 

data 

G4 

0.010+0.006 

0.018+0.001 

0.017+0.005 

0.040+.0007 

0.008+0.004 

0.008+.0004 

0.009+0.002 

0.018+.0005 


data 

0.006+0.005 

0.001+0.003 

0.28 


G4 

0.014+0.001 

0.007+.0003 


Data 

0.023+0.006 

0.017+0.003 

1.06 


G4 

0.065+0.001 

0.029+.0007 


0.004+0.003 

0.003+0.002 

0.0026+0.0026 

0.004+.0002 

0.002+.0001 

1.8E-05+1.8E-06 

0.006+0.002 

0.002+0.002 

0.0007+0.0007 

0.009+.0003 

0.004+0.002 

3.6E-05+1.8E-06 

0.004+0.004 

0.003+0.003 

0.00017+0.00017 

0.003+.0002 

0.002+.0001 

1.8E-05+1.8E-06 

0.007+0.003 

0.002+0.0031 

0.0024+0.0024 

0.013+.0004 

0.007+0.0003 

7.2E-05+0.2E-05 


CH2 0.99 

C 1.53 

A1 1.69 

Fe 0.82 


data 0.023±0.005 0.013±0.002 0.008±0.003 0.002±0.002 0.0005±0.0005 

G4 0.054±0.001 0.025±.0006 0.012±.0004 0.006±.0003 5.4E-05±1.8E-06 

data 0.032±0.008 0.015±0.003 0.007±0.002 0.004±0.001 O.OOlliO.OOll 

G4 0.061±0.001 0.026±.0006 0.013±.0004 0.007±.0003 5.4E-05±1.8E-06 

data 0.071±0.009 0.028±0.003 0.018±0.003 0.005±0.002 0.0008±0.0008 

G4 0.124±0.002 0.055±0.001 0.026±.0006 0.013±.0004 0.00010±3.6E-06 

data 0.059±0.008 0.026±0.003 0.015±0.003 0.006±0.002 0.0014±0.0013 

G4 0.099±0.001 0.045±.0008 0.022±.0005 0.011±.0004 9.0E-05±1.8E-06 



data 

0.022+0.006 

0.013+0.003 

1.31 

G4 

0.071+0.001 

0.036+.0008 


ELMS 

0.021+.0009 

0.010+.0006 


data 

0.039+0.009 

0.017+0.004 

1.90 

G4 

0.128+.0015 

0.052+.0009 


ELMS 

0.034+0.001 

0.016+.0007 


0.008±0.004 0.007±0.005 0.0025±0.0025 

0.023±.0006 0.017±.0004 0.0024±1.4E-05 

0.005±.0004 0.002±.0002 1.8E-05±1.8E-06 

0.010±0.002 0.008±0.003 0.004±0.004 

0.029±.0006 0.019±.0005 0.0025±1.5E-05 

0.007±.0004 0.003±.0003 3.6E-05±0.2E-05 


Table 2 

The data and simulated predictions for the probability per radian of scattering 
through a given projected angle. The table is split on two halves to fit the page. 
The first bin starts at zero. The errors on each bin are total errors. Note that the 
last bin is essentially an overflow bin and has 100% errors. 
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